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Abstract. - We calculate the mesoscopic fluctuations of the magnetic anisotropy of ferro- 
magnetic nanoparticles; that is, the effect of single-particle interference on the direction of the 
collective magnetic moment. A microscopic spin-orbit Hamiltonian considered as a perturbation 
of the much stronger exchange interaction first yields an explicit expression for the anisotropy 
tensor. Then, assuming a simple random matrix model for the spin-orbit coupling allows us 
to describe the fluctuation of such a tensor analytically. In the case of uniaxial anisotropy, we 
calculate the distribution of the anisotropy constant for a given number of electrons, and its 
variation upon increasing this number by one. The magnitude of the latter scales inversely 
with the number of atoms in the particle and is sufficient to account for the experimental data. 



In a remarkable step toward a better understanding of ferromagnetism at the nanometer 
scale, individual electronic states within ferromagnetic nanoparticles were recently observed 
using electron tunneling spectroscopy [1,2]. Several features of the tunneling excitation spec- 
tra could not be explained within the simple electron-in-a-box picture, unlike the normal metal 
case [3] . Although this "failure" was anticipated due to the many-body character of ferromag- 
netism, several particular aspects were largely unexpected: (1) the non-linear magnetic field 
dependence of the energy levels close to the switching field (B SV! , the field at which the magne- 
tization suddenly changes its orientation), (2) very small g- factors, and (3) excitation energies 
much smaller than the mean level spacing associated with the finite particle size. While we 
briefly return to the latter two in closing, we focus on the first effect as representative of the 
peculiarities of nanoscale ferromagnetism. 

A non-zero B sw implies the presence of a magnetic anisotropy energy barrier: the (large) 
magnetic moment of the particle points in a preferred direction. This preference originates 
from spin-orbit coupling and/or the shape of the grain. At B sw , there is an abrupt change in 
the grain's total energy, leading to "jumps" in the tunneling spectrum. Refs. [1,2,4] show that 
such jumps can be described by a simple uniaxial anisotropy model: Ha = —kS^/S, where 
k is the anisotropy constant — which is an intensive quantity in the bulk — S is the total spin 
and z is the anisotropy axis. However, to account for the non-linear behavior and the magni- 
tude of the jumps, it was crucial to assume that k depends on both the number of electrons 
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in the nanoparticle and the particular many-body states involved in each electronic transi- 
tion. A change of 1-3% upon adding one electron in a nanoparticle with roughly 1000 atoms 
was enough to account for the experiments. Although fluctuation of electronic properties is 
expected at the nanoscale, a microscopic connection to anisotropy is needed. 

In the last years, several works have aimed to understand the origin of the additional 
resonances [2,4-7] or the fluctuation of the anisotropy [8,9]. While the former remains con- 
troversial, the latter seems to be accounted for (in nanoparticles with up to 260 atoms) by the 
self-consistent numerical calculation done in Ref. [8]. One may wonder, however, if a simpler 
model, suitable for an analytic treatment of larger nanoparticles, can give similar results while 
providing additional physical insight. 

In this work, we first derive an expression for the magnetic anisotropy tensor in terms of 
the matrix elements of the spin-orbit (SO) Hamiltonian using a simple Hamiltonian model 
for the ferromagnetic nanoparticle. Then, by using random matrix theory (RMT) to account 
for the fluctuations of the SO matrix elements, we calculate the probability distribution of 
the magnetic anisotropy constant for the case of uniaxial anisotropy. For nanoparticles with 
N a ~ 1000 atoms, we find that fluctuations of 1% are typical. Furthermore, the fluctuations 
scale inversely with N a so that our results agree with those of Ref. [8] for smaller nanoparticles. 

Ferromagnetism originates in the electron-electron interaction — in particular, the exchange 
interaction. The simplest model Hamiltonian, then, reads 

H = H s +H z + H so =J2 £ m ™ m o ~ JsS 2 + H Z + H so , (1) 

m,cr 

where {e m } are the single particle energies, S is the total spin operator, and Js is the exchange 
constant. The latter produces the large spin of the nanoparticle (typically S ~ 1000 [1]); notice 
that this model assumes that the magnetic moment is coherent and uniform [8]. The second 
and third terms in Eq. Q describe the Zeeman energy, Hz = — g/XB-B-S 1 , and the SO coupling, 

HSO = i ^ nm ' daa ' ) C na c ma' > ( 2 ) 

respectively. Hso is responsible for the magnetic anisotropy — we neglect the contribution of 
the shape anisotropy (magnetostatic) as its relative contribution decreases for smaller particles 
and so is very small at the nanoscale [10]. Here, it is a vector of Pauli matrices, and 

Aura = -l-^^i^Ljlv) ll) n (j, H)ljj m (j,u) , (3) 

where j labels the atomic sites and /x, v the orbital bands (s, p, d). 4> n (j, jJb) — (j^\ip n ) is the 
n th single-particle wavefunction, and Lj is the angular momentum with respect to site j. The 
parameter £ characterizes the SO magnitude; it is of order 80 meV for cobalt [8]. The SO 
matrix elements satisfy A nm — —A mn due to time-reversal symmetry. We consider the case 
where both the exchange energy and the single particle level spacing are large compared with 
the SO matrix elements (true for small enough nanoparticles) and treat Hso as a perturbation. 
In addition, because cobalt nanoparticles apparently have a fee structure [10], we do not 
include here the crystalline magnetic anisotropy of bulk hep cobalt. 

We start by deriving an effective spin Hamiltonian for the nanoparticle. In the absence 
of SO coupling both the occupation numbers and the total spin are good quantum numbers. 
In its presence, only states with different occupation numbers are mixed (recall Hso is fully 
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Fig. 1 - The single particle energy levels of Hs for N and TV + 1 electrons. The arrows point out 
different classes of transitions induced by the SO coupling, labeled I, II, III, and IV. Here <5 m (8m) is 
the mean level spacing for the minority (majority) spins at the Fermi energy. SEf is the single-particle 
energy difference between the highest occupied levels of the majority and minority spins. 



off-diagonal, A nn = 0). Therefore, it is natural to write an effective Hamiltonian restricted 
to a subspace of fixed occupation which includes virtual transitions out of such subspace due 
to SO coupling (up to second order). This is done by a suitable canonical transformation, 
H c ff = cxp[— \U]H cxp[i?7] . Let {\q)} denote the eigenstates of H$, Hs\<l) = E q \q). We choose 
U such that U q / q = if \q) and \q') have different occupation numbers and 



U q > q 



Va> 



{E q - - E q ) 



(4) 



if they do not [11]. Here V q > q = (Hso)q'q and B q > q is proportional to Hz and Hso- As for the 
case of normal metal grains [12, 13], the latter leads to a correction to the g-factor which we 
absorb into the Zccman term. Then, to leading order 



{H c ff) q ' q = E q Sq'q 



(Hz) q > q + lj2 V ^" 



q"q 



E„ — E„ 



Eq' " E q n 



(5) 



if \q) and \q') have the same occupation numbers, and (H e s)q>q — if not. Using the fact 
that the exchange interaction is very large (5Ep ^> V q ' q ), Eq. ||SJ can be further restricted 
to a subspace of a given total spin S. After using the Wigner-Eckart theorem [14] and some 
tedious algebra, we finally obtain the following spin Hamiltonian 



H e s = Hs + H' z 



y / J S a K. a pSp/S. 



(6) 



with an explicit expression for the anisotropy tensor 7C given below. Note that while H c h 
and H have the same eigenenergies, they do not have the same eigenvectors — in fact, \E) = 
exp [i£7] | E g) . Therefore, while the eigenstates of Eq. © have well-defined spin and occupation 
numbers, this is not the case for the eigenstates of the physical Hamiltonian (|TJ, as expected 
due to SO coupling. The spin selection rules for electron tunneling are thus less restrictive 
than for a state with well-defined spin. 
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We classify the transitions induced by the SO coupling in four types (see Fig. 0: from 
doubly-occupied levels to singly-occupied (I) or empty (II) levels, and from singly-occupied 
to singly-occupied (III) or empty (IV) levels. With this notation, the anisotropy tensor reads 



-2J S (6 I +6 W ) 



(e n — e m ) {e n — e m — 2JsS) 



4J S S 



(En £m) 2 Aj"gS 2 



25 



in 



(25-1) 



12 1 2S 




[En-Em (S + l) E n - S m - 2 J S S {2S + 1) E n -E m + 2 J S (S + 1) (5 + 1) (25 + 1 



where Sx, X = I, II, III, IV, restricts the sum to the appropriate set of transitions. It is 
clear from this equation that states with different occupation numbers will have different 
anisotropy tensors, leading to level-to-level fluctuations of the anisotropy constant. We then 
conclude that the addition of an electron to the nanoparticle will produce a different anisotropy 
tensor depending on the new occupied level. 

So far we have not specified any features of the SO coupling: Eq. (0 is valid provided only 
that the SO coupling is small. A detailed treatment of the anisotropy tensor for a particular 
grain geometry would require the calculation of the SO matrix elements using, for instance, 
a non-interacting tight-binding model or a more sophisticated self-consistent Hartree-Fock 
approach [8]. However, our goal here is to calculate the generic features of the fluctuations 
of the anisotropy tensor. Since the shape of the grain is assumed to be irregular, we expect 
the SO matrix elements to fluctuate from level to level due to the fluctuations of the single- 
particle wavefunctions. The usual RMT approach of describing the single-particle fluctuations 
with the Gaussian unitary ensemble then implies that the {A nm } are Gaussian uncorrelated 
random variables because of the sum over a large number of terms in their definition Eq. J2J . 
This assumption of randomness of the SO matrix elements is consistent with the results of 
Ref. [8]. On the other hand, we know from Refs. [2,4] that a uniaxial anisotropy is needed to 
describe the main features of the non-linear behavior in the experiment. Therefore, we assume 
uniaxial anisotropy on average: (A 2 ) = (A 2 ) and (A 2 ) < (A 2 ). In a more detailed calculation, 
this asymmetry would probably come as a difference in the energy dependence of the SO matrix 
elements (this is the case in Ref. [8]); in our model, however, it can only appear as a difference 
in the standard deviations. The case of pure "mesoscopic anisotropy" — i.e. when all (A 2 ) are 
equal — was recently analyzed in Ref. [9]; it was found that for N a > 100 mesoscopics alone 
cannot explain the full magnitude of the anisotropy, thus indirectly supporting our insertion 
of an additional uniaxial anisotropy. Since we do not expect (A 2 ) to be much smaller than 
(A 2 ), we take (A 2 } = (A 2 }/2; our results are only weakly sensitive to this choice (see below). 

In the basis where 7C is diagonal, the anisotropy Hamiltonian takes the form 

Ha = ^S 2 + (k 3 - ^)5f/5 (8) 

where {fci, fe, fo} are the eigenvalues of 7C and 53 = 5- £3. The first term is a small correction 
to the exchange constant and can be neglected. The second term represents the uniaxial 
anisotropy. We have verified that for the values used below for the SO matrix elements, 
(z-k 3 ) ~0.99, so that Ha — — kjyS 2 /S with fcjv = {k% + fc 2 )/2 — fc 3 . The fluctuations of fcjv 
are calculated by numerically diagonalizing K, for each realization of {-A" m }. No fluctuation 
is introduced for the energy levels — they are taken to be locally equally spaced, in such a way 
that they reproduce the (hep) Co density of states. 

The distribution of kjy for a nanoparticle containing N a = 1200 atoms is shown in Fig. 
121 where we have used the following values for the parameters: SE F — 2eV, 6m = 3.7meV, 
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Fig. 2 - Probability distribution of the magnetic anisotropy constant for a nanoparticle size of N a = 
1200. The distribution is Gaussian with (/cat) = 23.7/^eV and cr(fcjv) = 0.4 /j,eV. 

5 m = 1.15meV, Js — 0.95meV, and a(A x ) = 0.32meV. The latter is estimated from Eq. 
J2Il using random and uncorrelated wavefunctions, (A 2 ,) — [£ 2 / '27 N a ](N s / N a ) — the additional 
factor N s /N a , with N s the number of atoms at the surface, is discussed below. The distribution 
is Gaussian with a mean value (k^) = 23.7 fj,eV and a standard deviation c(fcjv) = 0.4 fieV, 
ie. cr(fcjv)/(fcAr) — 2%. These values can be estimated from Eq. by replacing the sums by 
integrals and considering only the diagonal terms of 7C. We then get 



where a and b are numerical factors that depend on the shape of the density of states. As 
expected, the fluctuations go to zero in the bulk limit, er(A;jv)/(fc/v) oc 5 m /5Ep oc -A" 1 up to a 
logarithmic factor. We also note that in bulk (A 2 ) ocA" 1 so that (fcjv) is independent of the 
system's size. However, since we assume the anisotropy comes from the breaking of symmetry 
at the surface of the nanoparticle, we take (A%) oc N s /N 2 , so that (fcjy) oc N s /N a [10]. For 
Co, we have a ~ 0.8 and b ~ 1.6. Evaluation of Eqs. J^J and l|10|) gives (fcj\r) — 29 /xeV and 
a(k]\r) — 0A fieV, in good agreement with the numerical data. 

Let us now calculate the change of the anisotropy constant upon the addition of one 
electron with the minority spin (right panel in Fig.^). We calculate fcjv+i as before but take 
into account that both fcjv and k^+i contain contributions from the same SO matrix elements. 
We obtain a Gaussian distribution for Sk^ = fc_/v+i— ^n, with (Skjsr) — —0.03 /xeV and a{5k^) = 
0.2/ieV. Therefore, for a nanoparticle containing ^^ 1200 atoms, <j(5fcjv)/(fciv) ~ 1%- This 
is of the order needed to account for the experimental data of Ref. [2] . 

Figure 01 shows a(SkN)/(kjsr) as a function of N a . The results for the experimental data 
and the self-consistent calculation are indicated by the square and the diamonds, respectively. 
It is important to emphasize that once the ratio (A 2 Z ) / (A^.) is fixed, and so is (fcjv), there 
are no free parameters left to fit the fluctuations. Thus, the fact that the relative fluctuation 
of the anisotropy constant agrees reasonably well with both the experimental data and the 
results of a self-consistent numerical study is quite remarkable and indicates that our simple 
model is able to capture the essential physics. 
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Fig. 3 - Relative fluctuation of the anisotropy constant upon the addition of one electron as a function 
of the numbers of atoms in the nanoparticles. Fluctuations are of the order of 1% for N a ~ 1000. 
Also shown are results of experiment (square, size of error bar corresponds to nanoparticles of 3-4nm 
diameter) [2,15] and of a self-consistent calculation (diamonds) [8]. 



Following the same procedure as above, we estimate 



(Sk N )^c[(Al)~(Al)] 



Js 



S m SE F 



(11) 



a(6k 



2J S 
S m SE F 



2[{Al)* + 2{Alf 



li 



- 1-4 J S 



<>m - Js 



(12) 



where c contains logarithmic corrections in SE F /S m . In our case c: 
-0.07 neV and cr(<5fcAr)~0.3/ieV. Note that 



1 



-3.5, which gives {5k 



<j(Sk N ) S m 

(k N ) ^ SE F K N a 



(13) 



also decays to zero in the bulk limit. We notice that the ratio a{5kpj)/ (kpf) is independent of 
the value of (A%.) and only weakly dependent (up to a factor 2) on the ratio {A 2 Z )/{A^,) if the 
latter is < 0.5. The same is true for the relative fluctuation of the ground state anisotropy, 
c(fcjv) / '{kN) ■ Consequently, our estimate of the magnitude of the fluctuations is quite robust, 
being largely controlled by the ratio 6 m /6E F . 

It is straightforward to check from Eqs. i|ll|) and (|12f) that a(5k]y)^> {Sk^). This implies 
that fcjv+i can be larger or smaller than fcjv with roughly equal probability. While this is a key 
point for a correct account of the experimental data [2] , it also has an interesting consequence 
if the dwell time for electrons in the nanoparticle is larger than the time it takes for the 
magnetization to switch. In that case, it should be possible to switch the magnetic moment 
of the nanoparticle by first tuning the magnetic field just below B sw for the ground state and 
then increasing the bias voltage until a high energy state with a smaller anisotropy constant 
is excited. The reverse procedure should present hysteresis. Gated nanoparticles should show 
similar effects as a function of the gate voltage. 

Finally, let us briefly comment on the other two features of the experimental data men- 
tioned in the introduction. First, as noted above, B q / q gives a correction to the g-factor 
[cf. Eq. As for the magnetic anisotropy, this correction will depend on the particular 
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electronic state of the nanoparticlc. We note that even very small fluctuations will strongly 
modify the effective g-factor, 



with SgN = (jN —Qn-i and AS = Sn — Sjv-i- For instance, Sgw/gN — 10~ 4 , S ~ 1000 and 
AS = — 1/2 can lead to the small g e s observed experimentally (this requires 5gjy/AS < 0, 
which should be verified). Second, regarding the excitation spectrum, the experimental data 
[15] suggest that the low energy spacing is roughly independent of the grain's size [16]. The 
only quantity with such scaling is the anisotropy constant. Then, resonances should be related 
to the excitation of collectives modes (magnetization) [7] along an axis with a much larger 
anisotropy (hard axis) — this is crucial since the smallness of B sw implies that the anisotropy 
along the easy axis is also small. We believe such resonances are described by a more general 
form of the anisotropy, such as Ha = —kSl + k'S%. Work in this direction is in progress. 
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